Vesicle dynamics in elongation flow: Wrinkling instability and bud formation. 
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We present experimental results on the relaxation dynamics of vesicles subjected to a time- 
dependent elongation flow. We observed and characterized a new instability, which results in the 
formation of higher order modes of the vesicle shape (wrinkles), after a switch in the direction of 
the gradient of the velocity. This surprising generation of membrane wrinkles can be explained 
by the appearance of a negative surface tension during the vesicle deflation, due to compression 
in a sign-switching transient. Moreover, the formation of buds in the vesicle membrane has been 
observed in the vicinity of the dynamical transition point. 

PACS numbers: 87.16.Dg, 47.60.+i, 87.17.Jj 



Giant unilamellar vesicles (GUV's), formed by a closed 
phospholipid bilayer, appear to be a well-defined simpli- 
fied system for studying physical aspects of the dynam- 
ics of biological cells. Equilibrium mechanical properties 
of vesicles are relatively well understood. Vesicles ex- 
hibit a variety of equilibrium shapes that correspond to 
minima of the membrane Helfrich's energy The non- 
equilibrium dynamics of a vesicle subjected to an external 
flow received intensive attention in numerous theoretical 
2, y, 3, la, B, numerical @, @, 0, LLQ|; an d experimental 
Til [T2I [~L3l 0, Ql| [l6[ studies. In simple shear flows a 
vesicle exhibits several types of motion - tank treading, 
tumbling and trembling fl3l [~L4j (also called vacillating- 
breathing Q or swinging |l(|), depending on its location 
in the space of the system control parameters (viscosity 
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contrast, excess area, shear rate) 

There is, however, a lack of experimental observations 
in the regime of transient dynamics, when the system 
undergoes a non-equilibrium relaxation toward one of 
its dynamically stable states. In stationary shear flows 
only the lowest order modes (usually just the second or- 
der modes) characterize the various dynamical states of 
the vesicles observed, since the energy contribution from 
higher order excitation modes, if present, would be much 
higher. However, as pointed out recently [17], the flow 
may happen to impose a negative tension on the mem- 
brane, which leads to the growth of the higher order 
modes and to a shape instability. 

Let us consider the Helfrich free energy functional in a 
general form [1]: 
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where k is the bending rigidity of the membrane, h is 
the local curvature and a is the vesicle surface tension, 
which is the Lagrange multiplier corresponding to the 
surface area conservation constraint. If we consider for 
simplicity just a fiat membrane, which can be param- 
eterized by a height function u(x : y), then the expan- 
sion of the functional F in the Fourier space up to sec- 
ond order in u gives: F^ = | ^ k (nek 4 + ak 2 )\uk\ 2 • 



One notices that modes Uk with k < ^\cf\/k become 
unstable for a < 0, resulting in the generation of the 
higher-order modes becoming energetically more favor- 
able. One possibility to experimentally realize a nega- 
tive surface tension of the vesicle membrane is the use 
of a time dependent flow, where the sign of the velocity 
gradients undergoes a fast change, under which a vesi- 
cle becomes temporarily deflated. The simplest realiza- 
tion of this idea is a plane elongation (hyperbolic) flow: 
v x = ex, v y = -ey, v z 0. 

In this letter, we present the first experimental study of 
giant vesicle dynamics in such time-dependent, transient 
plane hyperbolic flow. We study the vesicle relaxation 
towards a new stationary state in two cases: from an 
equilibrium state when the elongation flow is suddenly 
turned on, e(t) = H{t)e$, where H(t) is the Heaviside 
step function, and when the elongation flow is suddenly 
reversed, e(t) = sign(t)eo, i.e. v xx changes from — eo to €q. 
T he station ary, stretched state is known to obey D sat = 
V/15A/327T, with (j> = {0;tt/2} @, i, [T3|. Here D = 
(L — B)/ (L + B), L and B are the large and small semi- 
axis of the elliptical approximation of the vesicle cross- 
section and (f) is the inclination angle with respect to the 
x axis. We assume that the membrane is impermeable 
for the time scale of the experiment, with the excess area 
A = S/R 2 — 47r, where S is the total surface area of the 
vesicle and R its effective radius, defined via the volume 
V = InR 3 . 

Measurements of the vesicle dynamics were conducted 
in the vicinity of the stagnation point (v x = v y = v z = 0) 
via epi-fluorescent or phase contrast microscopy. The 
flow was produced in a cross-slot micro-channel of 500 
/am wide and 320 jam in height manufactured in elas- 
tomer (PDMS) by soft lithography [H[. The details of 
the design and of the arrangement will be published else- 
where. Particle tracking velocimetry measurements of 
the flow field show that the deviation of the elongation 
rate (Ae) xy /e across the size of the observation window 
is < 5%, deviations of e in the z-direction on the scale 
of the vesicle were (Ae) z /e < 5%, and that the ratio of 
shear velocity gradient j z to e on the size of the vesicle 
did not exceed (Ae) z /e. Experiments were performed in 
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FIG. 1: Relaxation dynamics in elongation flow suddenly 
turned on. Snapshots of a vesicle with (a) A = 1, X — 2-2, 
A = 0.5, e = 0.1, (b) A = 1, X = 4 - 4 , A = 0.6, e = 0.9. 
Numbers on the snapshots are the non-dimensional times te, 
the scale bar is 10 fim. (c) D versus te with filled squares from 
data in (a), open squares from data in (b). Inset: statistical 
average of T* versus A. 



the range of velocity gradients e = 0.05 -=- 10. We define 
the dimensionless strain x — tVoutR 3 / where n 10 -12 
erg for DOPC [19[. The viscosity of the fluid inside the 
vesicle, 7^ n , can be different from the viscosity rj out of the 
surrounding fluid, and their ratio A was varied across the 
experiments. The lipid solutions consisted of 85% DOPC 
(Sigma) and 15% NBD-PC (fluorescent lipid, Molecular 
Probes) dissolved in 9:1 v/v chloroform- methanol solvent 
(1.8 mg total lipids/ml), or DOPC in the solvent (1.5 
mg/ml). The methods and conditions of the preparation 
of the vesicles for the experiments have been described 
previously [II [II Bit 

The first set of the experiments was performed sud- 
denly switching on the flow, starting after the vesicle had 
relaxed into an equilibrium shape. An initial growth of 
D, monotonic in time till saturation to D = D sat was 
observed, whereas <\> reached the stationary value either 
of or 7r/2 (see FigJT]). We characterize the process for 

D « D sat by the linear growth time T* = e [^j] ~ . It 
was found that T*(x) « const in the range of x = 1 15 
for a given A. The dependence of T* on the viscosity 
contrast is found to be linear in the range A = 0.1 7 
(see inset FigJTh). Averaging is done on the data sets of 
30 to 100 points for each A. These experimental results 
are consistent with the recent theoretical predictions for 
the same relaxation to a stationary state for x > 1- The 
theory too shows that the relaxation time scales linearly 
with e and with A jUELEzl- 

In the second set of experiments, in which the flow was 
suddenly switched from v xx = — eo to v xx = eo ("sud- 
denly" essentially means that the switching time is much 



smaller than x/eo), vesicles undergo a relaxation from 
one stretched stationary state (D = D sat \(j) = tt/2) to 
another one (D = D sat \(j) = 0). Surprisingly, we found 
that above some value of x > Xc the vesicles develop 
small wavelength perturbations in the shape, which we 
call wrinkles (see Fig|2ji,b,c). 

Quantitative evaluation of the higher order modes of 
the membrane shape was performed in the following way: 
the vesicle contour is fitted by an ellipse and the am- 
plitude A(0, eot) of the deviation from the elliptical fit 
is taken as function of the angle for every instant of 
time (see upper inset of Fig|3j). Examples of A(6, cot) are 
shown in Fig|2]i. The instantaneous Fourier transforms 
of the amplitudes with respect to define the dynamics 
of the spectrum Uk(t). The temporal evolution of these 
power spectra is shown in Fig|2^. A difference in the 
development of the higher k modes for different values 
of x can be clearly seen: during the transition from one 
steady state to another, more higher order modes are ex- 
cited for larger \. The evolution of D is shown in Fig|2f 
in order to precisely define the time interval where the 
transitional dynamics takes place. 

The dependence of the average in time power spectrum 
of the relaxing modes, = |i£/c| 2 , for different values of 
X, is shown in the lower inset of Figl3j We found that 
the spectra show a oc k~ A dependence for x l ess than 
some critical value Xc? that is a well-known spectral de- 
cay due to thermal noise [1]. For larger x? the spectra 
become rather flat at smaller /c, while the higher modes 
still comply with the equilibrium spectral decay. From 
the rather sharp transition from the flat to the k~ 4 spec- 
trum around some k = kthr, we can postulate that for 
X > Xc modes with k < k t hr are excited dynamically, 
while the higher modes are excited not by the flow but 
rather by thermal noise. We determine Xc as the thresh- 
old above which the modes with k > 3 are excited (k = 3 
is the first mode higher than elliptical). To this extent, 

we define k* = \J^2kL 3 k 2 Pk/Y^k=3 The restriction 
to k < 19 is dictated by the image resolution of the small- 
est vesicles in the experiments. The dependence of k* on 
X, averaged over w 200 data points, is shown in Fig|3l 
For x < Xc, k* remains constant, which means that the 
spectrum of all modes with k > 3 obeys the equilibrium 
distribution, k~ A . The growth of k* starts at x = Xc 
with Xc = 6.5 ± 0.8, which is identified as the onset of 
the wrinkling instability. The dependence of k* above the 
instability threshold can be fitted by ~ x 1 ^ 4 m a good 
agreement with numerical simulations in the same range 
of x [17]. All the experiments were done for A = 1, and 
no investigation on the influence of viscosity contrast on 
the wrinkling phenomena was done. 

As we pointed out in the introduction, the generation 
of higher order modes in the vesicle shape, i.e. wrinkles, 
observed during the relaxation dynamics should lead to 
a tremendous increase in the elastic energy, which is un- 
likely, for a vesicle with positive surface tension. On the 
other hand, a recent theory predicts that a sudden switch 
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FIG. 2: Wrinkling instability. Snapshots of vesicle dynamics in time- dependent elongation flow at A = 1, A « 1 and: (a) 
X = 8.1, (b) x — 81, (c) x — 323.5. The scale bar is 20 /im, numbers are te. Plots below the images are the data analysis 
for each of the cases above: (d) - amplitudes A{6,t) of higher harmonics versus and te (values in color), (e) |^ife| 2 (£) are the 
instantaneous Fourier spectra of the amplitudes A(0,t) at various te (values in color), (f) D(t) versus te. Columns ex a , ex^, 
ex c correspond to the data presented in row (a), (b), (c). 



in direction of the velocity gradient, e(t) = sign(t)eo^ is 
effectively equivalent to a negative surface tension and 
leads to the excitation of modes with k 2 < /3x/VA, 
where (3 is the numerical factor [17[. Then the theory 
gives Xc — 1-2 for k = 3 and A c± 0.6, which corresponds 
to our average (A^) over the data set in the transient 
region. The theoretical value of Xc is of the same order 



as the experimental one, and the difference can be at- 
tributed, first, to uncertainty in the value of n taken and 
to rather rough estimates based on an isotropic surface 
tension [17j . 

Another interesting phenomenon observed is the for- 
mation of buds. These could be seen intermittently in 
the experiments: sometimes the vesicle surface folds to 
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FIG. 3: k* versus %. The arrow defines the onset of 
excitation of the mode k = 3. Dashed line is a fit ~ X 1 ^ 4 
above the instability threshold. The upper inset illustrates 
the image analysis. Lower inset: averaged power spectrum 
for various x : squares - 2.6; open circles - 24; circles - 116; 
dashed lines show a oc k~ A dependence. 




FIG. 4: Formation of buds, (a) x = 6.8, A = 0.9; (b) x = 179, 
A = 0.4. The scale bar is 10 fim, numbers are te. 

the point that it creates the enclosure of a smaller vesicle 
inside the main one (FigHJ). We have not studied this 
phenomenon in detail, but some of its features can be 



described as follow: the enclosure process is irreversible, 
the bud scale is much smaller than the scale of the excited 
mode corresponding to the given x, and the phenomenon 
is mostly seen in the vicinity of Xc- 

To summarize, we presented new experimental results 
about the relaxation dynamics of vesicles in elongation 
flows suddenly switched on or reversed. When the vesicle 
relaxes from its equilibrium shape towards a new station- 
ary, stretched shape in the elongation flow, the scaling of 
the dynamics with the elongation rate and the linear de- 
pendence of the relaxation velocity on the viscosity con- 
trast were found to be in agreement with the theoretical 
predictions 0, @, [TtJ • We also observed and characterized 
a new instability, that results in the excitation of higher 
order modes, i.e. wrinkles, in the membrane, during the 
vesicle relaxation following the reversal of the velocity 
gradient. This unexpected generation of higher order 
modes suggests that only the appearance of a negative 
surface tension during the vesicle deflation due to com- 
pression in the transient can explain the effect. A recent 
theory [17[ used this physical picture to derive a crite- 
rion for the onset of the instability and the power law 
dependence of the average wave number of the higher 
order modes as a function of x> which is in reasonable 
agreement with our experiment. Finally we observed and 
report here, albeit without a quantitative investigation, 
the phenomenon of bud formation during the transient 
dynamics, particularly close to the wrinkling instability 
threshold. 
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